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Phase space factors and half-life predictions for Majoron emitting (5 jS decay 

J. Kotila,^’^’0 J. Barea,^’0 and F. lachello^’d 

^ University of Jyvaskyla, Department of Physics, 

B.O. Box 35, FI- 40014 , University of Jyvaskyla, Finland 
^Center for Theoretical Physics, Sloane Physics Laboratory, 

Yale University, New Haven, Connecticut 06520-8120, USA 
^Departamento de Fisica, Universidad de Concepcion, Casilla 160-C, Concepcion 4070386, Chile 

A complete calculation of phase space factors (PSF) for Majoron emitting QvP~I3~ decay modes 
is presented. The calculation makes use of exact Dirac wave functions with finite nuclear size and 
electron screening and includes life-times, single electron spectra, summed electron spectra, and 
angular electron correlations. Combining these results with recent interacting boson nuclear matrix 
elements (NME) we make half-life predictions for the the ordinary Majoron decay (spectral index 
n=l). Furthermore, comparing theoretical predictions with the obtained experimental lower bounds 
for this decay mode we are able to set limits on the effective Majoron-neutrino coupling constant 


PACS numbers: 23.40.He, 23.40.Bw, 14.60.Pq, 14.60.St 


I. INTRODUCTION 

Double-/? decay is a process in which a nucleus {A, Z) 
decays to a nucleus {A, Z ± 2) by emitting two electrons 
or positrons and, usually, other light particles 

{A, Z) — {A, Z ± 2) -I- 2e^ -I- anything. (1) 



The mode where two antineutrinos or neutrinos are emit¬ 
ted is predicted by the standard model and has been ob¬ 
served in several nuclei (for a review, see e.g. Q)- The 
more exotic mode, neutrinoless double beta decay, is not 
allowed by the standard model, and once observed would 
offer new information on many fundamental aspects of 
elementary particle physics. As discussed in Ref. [H, sev¬ 
eral scenarios of neutrinoless double beta decay have been 
considered, most notably, light neutrino exchange, heavy 
neutrino exchange, and Majoron emission. After the dis¬ 
covery of neutrino oscillations, attention has been focused 
on the first scenario and the mass mode, where the tran¬ 
sition operator is proportional to {mj) jraf,. Even though 
most current experimental efforts have been focused to 
the detection of this mode, interest on the mechanism 
predicting 0;z/3/3 decays through the emission of addi¬ 
tional bosons called Majorons has also renewed lately. 

Majorons were introduced years ago [1,0| as massless 
Nambu-Goldstone bosons arising from a global B — L 
(baryon number minus lepton number) symmetry broken 
spontaneously in the low-energy regime. These bosons 
couple to the Majorana neutrinos and give rise to neutri¬ 
noless double beta decay, accompanied by Majoron emis¬ 
sion OuPffM as shown in Fig. [T] (a). Although these 
older models are disfavored by precise measurements of 
the width of the Z boson decay to invisible channels Q, 
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FIG. 1. Neutrinoless double beta decay accompanied by the 
emission of one or two Majorons. 

several other models of QvfffiM decay have been pro¬ 
posed in which one or two Majorons, denoted by xo, are 
emitted, (see Fig. [1]) 

(A, Z) —(A, Z -h 2) -h 2e + xo (2) 

or 

(A, Z) —> (A, Z -h 2) -h 2e -I- 2xo- (3) 

Table U lists some of the models proposed to describe 
these decays. The different models are distinguished by 
the nature of the emitted Majoron(s), i.e. is it a Nambu- 
Goldstone boson or not (NG), the leptonic charge of the 
emitted Majoron(s) (L), and the spectral index of the 
model (n), which characterizes the shape of the summed 
electron spectrum, as described in Sect. El 


In our previous articles we have studied phase space 
factors (PSF) and prefactors (PF l fl2-ll4. and nu¬ 
clear matrix elements (NME) [2, Il3l. fl5l - ll^ needed 
for the theoretical description of 2^/3“/3“, 

0J//3+/3+, OvP+EC+, RQvECEC, 2i//3+^+, 2vl3+EC, 
and 2vECEC decay mediated in the case of neutrino¬ 
less decay by light or heavy neutrino exchange. In this 
article we continue our systematic evaluation by present¬ 
ing phase space factors for the different Majoron emit- 
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TABLE I. Different Majoron emitting models of QvPfiM de¬ 
cay IM- The third, fourth, and fifth columns indicate 
whether the Majoron is Nambu-Goldstone boson or not, its 
leptonic charge L, and the model’s spectral index n. 


Model 

Decay Mode 

NG boson 

L 

n 

IB 

Oi4/3/3xo 

No 

0 

1 

IC 

QvPPxq 

Yes 

0 

1 

ID 

014/3/3x0X0 

No 

0 

3 

IE 

0y/3/3xoXo 

Yes 

0 

3 

IIB 

OvP/3xo 

No 

-2 

1 

lie 

OvP/3xo 

Yes 

-2 

3 

IID 

014/3/3x0X0 

No 

-1 

3 

HE 

0i4/3/3xoXo 

Yes 

-1 

7 

HE 

Ov/3/3xo 

Gauge boson 

-2 

3 

"Bulk" 

Ov/3/3xo 

Bulk field 

0 
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ting mechanisms and combining the n = 1 results with 
recent interacting boson model (IBM-2) NMEs [l^ to 
make half-life predictions for ordinary Majoron decay 
(n = 1). Furthermore, we compare our theoretical pre¬ 
dictions with the obtained experimental lower bounds 
for this decay mode to set some limits on the effective 
Majoron-neutrino coupling constant (g^). 


II. PHASE SPACE FACTORS IN MAJORON 
EMITTING DOUBLE-d DECAY 


The key ingredients for the evaluation of phase space 
factors in single- and double-/? decay are the scattering 
wave functions and for EC the bound state wave func¬ 
tions. The general theory of relativistic electrons and 
positrons can be found e.g., in the book of Rose M- The 
electron scattering wave functions of interest in f5~ 
were given in Eq. (8) of [a. 

In order to calculate PSFs for Majoron emitting 
we use the formulation of Doi, Kotani, and Takasu 
The differential rate for the decay is given by (2^, 

dWmxon = C0S6»12) Wmxondeide 2 d{cos 0 i 2 ) 

(4) 

where ei and £2 are the electron energies, 0 i 2 the an¬ 
gle between the two emitted electrons, and Wmxon takes 
different values depending on the number of emitted Ma- 



jorons m and the spectral index n: 


Wlxol = 

WlXo3 = 
«^2xo3 = 

^2x0 r = 


g\{Gcos0c)‘^ f he 


QAir'^h 

g\{G cos Oc)'^ 3 


q{pic){p2c)eie2 


QAir'^h 

g\{G cos Oc)'^ 
30727r9h 

g\{G cos 0c)'^ 


\2R 
q^{pic){p2c)eie2 
f he'' ^ 


(WeC^) 


,2^-2 


V2i? 


g^(pic)(p2c)eie2 


(meC ) 


,2^-6 


he \ ^ 

:) q\pic){p2c)eie2. 


2RJ 


(5) 


Here G is the Fermi constant, 0c is the Cabibbo angle, 
and the Majoron energy q is determined as q = Qpp + 
2meC^ — Cl — £ 2 ) and R = with rg = 1.2 fm, is the 

nuclear radius.The quantities and in Eq. (|3]) can 
be written as (2lj 




M, 


(m,n) 


Oi^M 


* = 0 , 1 , ( 6 ) 


where | {g^) | is the effective coupling constant of the Ma¬ 
joron to the neutrino, m = 1,2 for the emission of one 
or two Majorons, respectively, and is the nuclear 

matrix element. The functions are defined as 

/!? = ir‘-‘l‘ + i/nl‘ + ir‘.t + l/r‘it , 

with 

= g-i{ei)g-i{e2), 

/ii =/i(ei)/i(e2), 

/-\=3_i(£i)/i(£2), 

/i ^ =/i(ei)ff-i(e2), 


where g-i{e) and fi{e) are obtained from the electron 
wave functions as explained in Ref. [l^ . 


All quantities of interest are then given by integration 
of Eq. ([1]). Introducing 

9 ^Q^/5-t-2meC^—ei 

( 7 O) =_ / / fO) 

g\\n2j^^,. (9) 

X Wmxondeide2j 

where the axial vector coupling constant gA is separated 
from the phase space factors for conveniency, we can cal¬ 
culate: 

(i) The half-life 



9a 


4 q (0) 


l(5x")l 




(m,n) 


Oi/M 


( 10 ) 
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(ii) the single electron spectrum 


dW, 


mxon _ iLr{m,n 


)dG. 


( 0 ) 

mxon 


de 


dei 


( 11 ) 


where ■ 

(iii) The summed electron spectrum, which shape makes 
the different Majoron emitting modes experimentally rec¬ 
ognizable 


d^^rnxon 

d{ei + € 2 ) 


Kf(n,in) 


J^(0) 

mxon 

d{ei + 62 ) ’ 


( 12 ) 


(iv) and the angular correlation between the two elec¬ 
trons 


a(ei) 


dG 

dG 


( 1 ) 

mxon 


( 0 ) 

mxon 


/dei 

/dei 


(13) 


TABLE II. Phase space factors Gm]cQn obtained using 
screened exact finite size Coulomb wave functions. 


Nucleus 


/^(O) /in 

^mxon\^^ 

1 

GO 

1 


m=I,n=I 

m=l,n=3 

m=2,n=3 

m=2,n=7 


1540 

17.1 

73.6 

690 

76Ge 

44.2 

0.073 

0.22 

0.420 

82Se 

361 

1.22 

3.54 

26.9 

96 zr 

905 

4.21 

11.0 

128. 

looMo 

598 

2.42 

6.15 

50.8 

llOpd 

94.1 

0.205 

0.487 

0.946 

ii^Cd 

569 

2.28 

5.23 

33.9 

i24Sn 

209 

0.653 

1.45 

4.45 

128Tg 

3.06 

0.001 

0.003 

0.0003 

ISOpg 

413 

1.51 

3.21 

14.4 

i34Xe 

2.92 

0.002 

0.003 

0.0002 

i36Xe 

409 

1.47 

3.05 

12.5 

148Nd 

197 

0.505 

0.986 

1.72 

isoNd 

3100 

21.1 

40.8 

538 

^®^Sm 

28.2 

0.034 

0.064 

0.021 

leoQd 

1590 

0.361 

0.672 

0.899 

198pt 

60.7 

0.068 

0.110 

0.021 

232Th 

82.4 

0.073 

0.105 

0.009 

238^ 

337 

0.532 

0.756 

0.213 


We have done a calculation of Gmxon and 
the list of nuclei shown in Table im We also plot our 
results in Figs. [2lf5l where they are compared with previ¬ 
ous calculations [201 - 1^ . For the comparison the values 
of [ 23 , [13 have been multiplied by a missing factor of 
two and divided by g\, and values of [2l| have been di¬ 
vided with factor g\AR^. The factor of two is the correct 
choice, as was acknowledged in [25|. This factor of two 
is also included in Table II of Ref. [13, where the de¬ 
cay of ^^®Xe have been studied. In their calculation for 


the phase space factors they use Fermi functions F(Z,E) 
that fully include nuclear finite size and electron screen¬ 
ing and are evaluated at the nuclear radius. Their result 
with different m = 1, 2 and n = 1,3, 7 are within 6 % of 
the ones reported here. 



40 60 80 100 120 140 160 180 200 220 240 


Mass Number 

FIG. 2. (Color online) Phase space factors in units 

(IQ-is yr“i). The label "approximate" refers to the results 
obtained by the use of approximate electron wave functions 
[20l - [^ . The figure is in semilogarithmic scale. 



FIG. 3. (Color online) Phase space factors in units 

(10“^^ yr~^)- The label "approximate" refers to the results 
obtained by the use of approximate electron wave functions 
[ 2 ^, [ 23 . The figure is in semilogarithmic scale. 

We also have available upon request the single and 
summed electron electron spectra and angular correla¬ 
tion for all nuclei in Table El An example, ^^®Xe decay, 
is shown in Fig.|51 The shape of the spectra is determined 
by the spectral index n, i.e. for the case n = 3 it is the 
same for both the emission of one or two Majorons, but 
the overall scale is different (yo result should be multi¬ 
plied by (/ic)^(2i?meC^)“^(487r^)“^ to obtain 2yo result). 
Also, since the angular correlation is obtained by divid¬ 
ing dGm\gn/dci by dGmxon/dci the calculation becomes 
unstable for the higher n near the endpoint energy and 
that region is thus excluded from Fig. [HI The distinc¬ 
tion between different values of n is most prominent in 
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FIG. 4. (Color online) Same as Fig. |3]for 
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FIG. 5. (Color online) Same as Fig. |3] for Due to 

a very small values, ^^®Te and ^^'^Xe are excluded from this 
figure. 


summed electron spectra, where as n increases the peak 
of the spectrum shifts from near the maximum kinetic en¬ 
ergy to near the minimum kinetic energy. The difference 
between different values of n is also shown in single elec¬ 
tron spectra but not as strongly as in the case of summed 
electron spectra. 


III. EXPECTED HALF-LIVES AND LIMITS ON 
COUPLING CONSTANT: ORDINARY 
MAJORON EMITTING d d" DECAY 

In the case of ordinary Majoron emitting double-/? de¬ 
cay, n = 1, the nuclear matrix elements have the same 
form as in the Qv(3l3 mediated by light neutrino exchange 

KIm = 9i[m^gt - (^-£) (14) 

The calculation of phase space factors can now be com¬ 
bined with updated nuclear matrix elements in IBM- 
2 [l^ to produce predictions for half-lives for ordinary 
Majoron decay {n = 1) in Table IIIII (left) and Fig. [3 


Judging by the predicted half-lives, the most prominent 
candidates are ^®°Nd and ^°°Mo. Furthermore, we can 
compare our predictions to half-life limits coming from 
experiments to set some limits on the effective Majoron- 
neutrino coupling constant. The obtained limits are 
shown on Table lllll (right). The most stringent limits 
for ordinary Majoron decay (n = 1) at the moment are 
for ^^®Xe coming from KamLAND-Zen [s^ and EXO- 
200 [ 2 ^ experiments reaching the order of magnitude of 
10-6. 


TABLE III. Left: Galculated half-lives for Majoron decay 
models where n = 1 (IB, IC, IIB) with = 10“"^, 

qa = 1.269 and recent IBM-2 nuclear matrix elements fi3| . 
Right: Upper limit on Majoron coupling constant (p^) from 
current experimental limits. 


8.19 

76Ge^76Se 39.8 

®2Se^82Kr 7.68 

96Zr^.96Mo 5.32 

looMo^iooRu 3.62 

iiopd_^iiocd 25.0 

ii6cd^ii6gn 7.06 

i28Te-)>i28xe 765 

i30Te_^i30xe 6.82 

i34Xe^i64gg^ gQ5 

i36Xe^i66Bg^ 20.1 

148Nd-^148Sjj^ 36.8 

i50Nd-).i50sm 1.74 

i64Sm^i64Qjj 273 

i60Gd^i60j3y 24.5 

198pt_^198Hg 232 

232Th_^232u 28.7 

238u_^238p^ 4.94 


> 7.2 X lO^o^ < 3.4 X 10-4 

> 6.4 X lO^^ < 7.9 X 10-6 

> 1.5 X 10^23 < 7.2 X 10-6 

> 1.9 X 1023 < 1.7 X 10-4 

> 3.9 X 10223 < 3.0 X 10-6 

> 8 X 1021^ 9.4 X 10-6 

> 2 X 1024_^ < 6.2 X 10-6 

> 1.6 X 102% < 6.5 X 10-6 

> 2.6 X 1024_^ < 6.2 X 10-6 

> 1.2 X 10243 < 9.2 X 10-6 

> 1.5 X 1021'= < 1.1 X 10-4 


Decay 


(102iyr) 


{g^e) (eV) 


48Ca^48Tj 



geochemical. 


IV. CONCLUSIONS 

In this article, we have reported a complete calcula¬ 
tion of phase space factors for Oi^/?-/?- decay proceeding 
through emission of one or two Majorons. The reported 
results include half-lives, single electron spectra, summed 
electron spectra, and electron angular correlations, to be 
used in connection with the calculation of nuclear matrix 
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FIG. 6. (Color online)Single electron spectra (panel a), summed electron spectra (panel b), and angular correlations between 
the two outgoing electrons (panel c) for the ^^®Xe —decay with different spectral indices n = 1, 3, 7. The scale 
of the panels (a) and (b) should be multiplied by for a realistic estimate. 


elements. Furthermore, we have combined our results 
with recent IBM-2 nuclear matrix elements to produce 
predictions of half-lives in the case of ordinary Majoron 
decay, spectral index n = 1. Comparing these predic¬ 


tions with experimental lower bounds, we have set some 
limits on the effective Majoron-neutrino coupling con¬ 
stant {g^). At the moment the best limits are com¬ 
ing from ^^®Xe experiments reaching the order of mag- 
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FIG. 7. (Color online) Expected half-lives for Majoron decay 
models where n = 1 (IB,IC,IIB) with = 10”“^, qa = 

1.269 and IBM-2 isospin restored nuclear matrix elements. 
The figure is in semilogarithmic scale. 


nitude of 10 Also, the results in Table III are for 
gA = 1.269. If Qa is renormalized to gA.eff, '’’ 1 / 2 ^ should 
be multiplied by (1.269/(7 a e//)"^ and limits on {g^) by 
{1.2m/gA.eff?. 
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